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Abstract — This paper investigates the sum-rate gains brought 
hy power allocation strategies in multicell massive multiple- 
input multiple-output systems, assuming time-division duplex 
transmission. For both uplink and downlink, we derive tractable 
expressions for the achievable rate with zero-forcing receivers 
and precoders respectively. To avoid high complexity joint opti¬ 
mization across the network, we propose a scheduling mechanism 
for power allocation, where in a single time slot, only cells that 
do not interfere with each other adjust their transmit powers. 
Based on this, corresponding transmit power allocation strategies 
are derived, aimed at maximizing the sum rate per-cell. These 
schemes are shown to bring considerable gains over equal power 
allocation for practical antenna configurations (e.g., up to a 
few hundred). However, with fixed number of users (N), these 
gains diminish as M —> oo, and equal power allocation becomes 
optimal. A different conclusion is drawn for the case where both 
M and N grow large together, in which case: (1) improved 
rates are achieved as M grows with fixed M/N ratio, and 
(ii) the relative gains over the equal power allocation diminish 
as M/N grows. Moreover, we also provide applicable values 
of M/N under an acceptable power allocation gain threshold, 
which can be used as to determine when the proposed power 
allocation schemes yield appreciable gains, and when they do 
not. From the network point of view, the proposed scheduling 
approach can achieve almost the same performance as the joint 
power allocation after one scheduling round, with much reduced 
complexity. 

Index Terms —Multicell, massive MIMO, power allocation. 


I. Introduction 

Multiple-input multiple-output (MIMO) systems are now 
regarded as a key technology for cellular communication 
systems due to various advantages they offer over single¬ 
antenna systems, including transmit diversity, high data rates 
and reliability [1-3]. By transmitting parallel data streams, it is 

Manuscript received November 15, 2014; revised March 30, 2015. 
This work was partly supported by the China 973 project under Grant 
2013CB329005, the China 863 Program under Grant 2014AA01A705 and the 
Jiangsu Graduate Research and Innovation Project under Grant KYLX_0806. 
The work of S. Jin was supported by the National Natural Science Foundation 
of China under Grant 61222102 and the International Science Sz Technology 
Cooperation Program of China under Grant 2014DFT10300. The work of M. 
McKay and D. Morales-Jimenez was supported by the Hong Kong Research 
Grants Council under grant number 616713. 

Q. Zhang and H. Zhu ai‘e with Jiangsu Key Laboratory of Wireless Com¬ 
munications, Nanjing University of Posts and Telecommunications, Nanjing, 
210003, P. R. China (email: zhangqiqi_1212@126.com; zhuhb@njupt.edu.cn). 

S. Jin is with the National Mobile Communications Research Labo¬ 
ratory, Southeast University, Nanjing, 210096, P. R. China (email: jin- 
shi@seu.edu.cn). S. Jin is the corresponding author. 

M. McKay and D. Morales-Jimenez are with the Department of Elec¬ 
tronic and Computer Engineering, Hong Kong University of Science and 
Technology, Clear Water Bay, Kowloon, Hong Kong (email:eemckay@ust.hk; 
eedmorales @ ust.hk). 


known that the ergodic capacity of MIMO channels increases 
linearly with the minimum number of transmit and receive 
antennas [4]. The transmitted information streams are sepa¬ 
rated at the receiver by means of appropriate signal processing 
techniques. The optimum receivers and precoders lead to a 
complexity burden on the system implementation. Therefore, 
linear receivers and precoders such as those based on the 
zero-forcing (ZF) criteria are often considered, as they offer 
significantly lower complexity with tolerable performance [5]. 

Recently, multiuser MIMO (MU-MIMO) systems, in which 
the antenna array of the base station (BS) simultaneously 
serves a multiplicity of autonomous users in the same time- 
frequency resource, have attracted substantial interest. Com¬ 
pared with conventional MIMO systems, MU-MIMO can offer 
a spatial multiplexing gain even if the users have only a single 
antenna [6-10]. 

Though MU-MIMO has been supported in 3GPP standards, 
it can not achieve the transmission rates demanded in 5G. 
Therefore, the concept of massive MIMO, which considers 
the use of hundreds of antenna elements to serve tens of 
users simultaneously, has come to the forefront of wireless 
communications research [11]. Massive MIMO can reap all 
the benefits of conventional MIMO at a greater scale, whilst 
providing substantial improvements in energy and spectral 
efficiencies with low complexity [12-23]. In particular, if 
the number of BS antennas is far greater than the number 
of single-antenna users, the simplest linear receivers and 
precoders become optimal [14]. It is also revealed in [13] that, 
when the number of BS antennas grows without bound, the 
achievable rate can be improved to a considerable level and 
the uncorrelated noise, fast fading and intracell interference 
all vanish. Moreover, [15] demonstrates that the transmitted 
power in massive MIMO systems can be proportionally scaled 
down with the number of BS antennas while maintaining the 
same rate performance. 

In power-limited systems, the equal power allocation policy 
cannot take full advantage of the energy resource and it is 
far from optimal. In this paper, allowing for different transmit 
powers at each user, we derive three tractable expressions for 
the achievable uplink rate—a lower bound, an upper bound 
and an approximation which lies between these two bounds. 
To avoid the joint optimization across the network, which 
involves high complexity algorithms with (potential) conver¬ 
gence issues, we proposed a scheduling mechanism for power 
allocation, where in a single time slot, only cells that do not 
interfere with each other adjust their transmit powers. Based 
on this, we subsequently derive power allocation strategies 


IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. X, NO. XX, XX 201X 


2 


which optimize each of these expressions and compare their 
relative performance. In addition, we present a tractable lower 
bound for the achievable downlink rate which is shown to 
be very tight, and similarly, a corresponding power allocation 
policy is derived. These schemes are shown to yield noticeable 
rate improvements over equal power allocation under various 
conditions. This is true, for example, for BSs housing up to a 
few hundred antennas (M). As M increases, the gain brought 
by optimized power allocation becomes less significant, with 
the optimal assignment approaching uniform allocation as 
M oo for fixed numbers of users (N). To further character¬ 
ize these effects more precisely, we also consider the regime 
where M and N grow large together, i.e., with M/N fixed, 
and show that; (i) improved rates are achieved as M and N 
both grow, and (ii) the relative gains over the equal power 
allocation depend on the ratio M/N, with larger gains for 
smaller ratios. By specifying a minimum required gain over the 
equal power assignment, we then provide the range of M/N 
ratios which satisfy the proposed constraint, and these results 
can be used to provide engineering insight into when the gains 
brought by optimized power allocation are sufficient to warrant 
the additional complexity relative to equal-power allocation, 
and when they are not. For example, for the uplink, if the total 
transmit power per cell is 20dB, systems with M/N < 12 (i.e., 
less than 12 BS antennas per user) can obtain more than 10% 
rate gain with optimized power allocation. These observations 
indicate that while equal power allocation is asymptotically 
optimal in the theoretical massive MIMO regime (i.e., when 
M/N —>• oo), for a wide range of practical scenarios in which 
the BS antennas is moderately larger than the number of users 
(but not infinitely larger), optimized power allocation can bring 
appreciable gains. 

This paper is organized as follows. Section HIl describes the 
multicell MIMO system model for both uplink and downlink, 
and dehnes the associated achievable rate with ZF receivers 
and ZF precoders, respectively. In Section |III] we present the 
lower and upper bound for the achievable uplink rate, followed 
by a more accurate approximation which lies between these 
bounds. Based on this approximation, the corresponding power 
allocation scheme is put forth. In Section |IV] we analyze the 
achievable downlink rate and, based on a tractable lower bound 
which we derive, an effective power allocation scheme is 
proposed. In Section |V] we provide a set of numerical results, 
while Section |VT] summarizes the main results of this paper. 

Notation —Throughout the paper, vectors are expressed in 
lowercase boldface letters while matrices are denoted by 
uppercase boldface letters. We use X^,X^,X* and X~^ 
to denote the conjugate-transpose, transpose, conjugate and 
inverse of X, respectively. Moreover, I^v denotes an iV x iV 
identity matrix, and [XJ^^ is the (i, j)th entry of X. Finally, 
E{-} is the expectation operator, || • || is the Euclidean norm, 
(a)+ denotes max{a,0} and tr (•) is the trace operation. 

II. System model 

Consider a multicell scenario with a central cell, denoted as 
cell 0, and L interference cells surrounding it, denoted as cell 
1 to L. Each cell contains a MU-MIMO system with N single¬ 


antenna users and one M-antenna BS. Channel reciprocity is 
exploited between uplink and downlink. 


A. Uplink 

On the uplink, it is assumed that all users simultaneously 
transmit data streams to their BSs in the same time-frequency 
resource and the transmit power for each user may be different. 
Let Gil {i,l = 0,1,..., L) denote the M x N MIMO channel 
matrix between the N users in the Zth cell and the M BS 
antennas in the ith BS. Therefore, the M x 1 received vector 
at the BS of the central cell can be written as 

L 

= Goof^y^xo + ^ Goi^y^xi + (1) 

1=1 


where the N x N diagonal matrix fti contains 
[P'1],... ,p^i,... ,p'^i] along its main diagonal, while 
p'//i is the transmitted power of the nth user in the hh cell, 
X; S denotes the transmitted signal vector from all 

users in the Zth cell, and iIq^ S represents the vector 

of additive white Gaussian noise with entries having zero 
mean and unit variance. 

We denote the independent channel coefficient between the 
nth user in the Zth cell and the mth antenna of the ith BS as 
gmini = [Gii]mn, which accounts for independent fast fading, 
geometric attenuation and log-normal shadow fading and can 
be expressed as [13] 


Sminl — ^minl \/l^inl, 


( 2 ) 


where hmini is the fast fading element from the nth user in 
the Zth cell to the mth antenna of the ith BS, which has 
independent real and imaginary parts with zero mean and 
variance 1/2. The large scale fading from the nth user in 
the hh cell to the ith BS models both the geometric attenuation 
and shadow fading. It is reasonable to assume the large scale 
fading coefficient to be constant across the antenna array, since 
the distance between users and the BS is much larger than the 
distance between antennas, and the value of fiini changes very 
slowly with time. Then, 

Goi = Ho^Do/^/^ (3) 


where Hq; denotes the M x N fast fading matrix between 
the users in the /th cell and the central cell’s BS, i.e., 
[Ho/]mn = hmOni and Doz is a iV X iV diagonal matrix 
with [Dozjnn = Poni- We assume all BSs have perfect 
channel state information (CSI) of channels in their own cell, 
which is a reasonable approximation in environments with 
low or moderate mobility. Linear processing is assumed at 
all receivers. Let Aqo be the M x N linear receiver matrix 
of the central cell’s BS, which depends on the channel matrix 
Goo- The received signal vector at the BS is processed as 

ro = (4) 

Lor ZL receivers, Aqo = Gqo (G^Gqo) • Substituting this 
along with ([T]i into (IDi gives 

L 

ro = -f ^ A^Gozrjy^xz -f A^rio^ (5) 
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The nth element of Tq can be further expressed as 


rnO = 




^noSOclXcl + a 


H 

OnO 


^0 ’ 


1^1 C=1 


(6) 


where ao„o is the nth column of Aqo, god is the cth column 
of Go;, and Xd denotes the cth element of the signal vector 
X; (c = 1,..,, iV; / = 1,,.., L), The desired signal in 
I© is y/p^Xno, and the remaining terms are considered as 
interference and noise. The interference-plus-noise term is 
modeled as Gaussian noise, which constitutes a “worst case” 
noise assumption [24], The ergodic achievable uplink rate of 
the nth user in the central cell is then given as follows [15] 


RnO = E 


log: 


1 + 


Ko 


L N 9 0 

E Pd\i^noSoci\ +||ao„o|| 


\ ;=ic=i 


J) 


(7) 


B. Downlink 

On the downlink, we assume that the transmit power sent 
by the BS for each user may be different. Then, the received 
signal by users in the central cell can be given as 

L 

= G^,Boo^l^\o +Y. + n^\ ( 8 ) 

1=1 

where B;; {I = 1,.,.,L) is the M x N precoding matrix 
applied at the /th BS, which is dependent on the channel 
matrix G;;, is a N X N diagonal power matrix containing 
,Pni] along its main diagonal, while jg 
the transmitted power of the Zth BS intended for the nth 
user. Moreover, q; € is the vector of signals intended 

for the N users in the Zth cell and iIq^ contains white 
complex Gaussian noise with entries having zero mean and 
unit variance. To satisfy the power constraint at the BS, B;; 
is chosen such that [25] 

E{tr(B„B/f)} =1. (9) 

We assume users have perfect CSI and detect received 

signals with the optimal maximum likelihood receivers. The 
ZF linear precoder is given by 

B„ = a;Grz(G^Gr,)“\ (10) 

where the constant scalar ai is chosen to conform to the 
constraint (|9]l. The substitution of (ITOl i into ® yields 

L 

= ao’S'o^'qo T^GfoB^vl/V^q, + (H) 


1=1 

and the nth element of yg^ can be further expressed as 

L N 



1=1 c=l 


where g;„o bi^ represent the nth and cth column of 
matrix Gjg and B^, respectively, qi^ is the cth element of the 
signal vector q; and n^g is the nth element of Gaussian noise 


vector ng^. Similarly with the uplink, the ergodic achievable 
downlink rate of the nth user in the central cell can be 
expressed as 



/ 

\ 

2 dl 

l0g2 

1 + 

aoPnO 

L N „ 


1 

T, Pd \sTn0^lcl\ +1 

c=l / 


. (13) 


C. Per-cell Optimization Approach 

Based on the achievable rates obtained above, we wish to 
design a power allocation scheme to optimize the network 
performance, A global (joint) optimization over all cells would 
be desirable, but it is found to be both very challenging and 
not practical due to; 1) high complexity of associated iterative 
algorithms (based on, e,g,, game theoretic ideas), which may 
have convergence issues as well, and 2) the unbounded and 
continuous (in space) deployment of BSs which makes it im¬ 
possible to define an isolated cluster (set of BSs) without being 
subject to other (out-of-cluster) interference. Therefore, we 
adopt a more practical but sub-optimal per-cell optimization 
approach which focuses on a single-cell performance while 
signals from other cells are regarded as constant interference. 
More precisely, the sum-rate of a (target) cell is optimized 
while the interference footprint (caused by surrounding cells) 
remains unaltered. This approach entails a certain (somewhat 
simple) level of coordination among cells, according to which 
power allocation is only performed by a cell when being 
granted (scheduled) to do so at a particular time slot. A 
network management entity may, for instance, grant different 
time slots to different cells according to some scheduling 
scheme, as exemplified next. 

Consider power allocation being applied in an scheduled 
manner to a set of cells. Cells that do not interfere with each 
other can perform power allocation simultaneously. Therefore, 
the entire set of cells can be divided into several groups, 
where each group consists of geographic distinct cells that 
do not interfere with each other and can, thus, complete 
power allocation simultaneously (i.e., similar to the adoption 
of frequency reuse patterns). Different groups are assigned 
with different time slots to allocate their transmit powers. 
Thus, within one scheduling round, each cell has operated 
with maximized sum rate for at least one time slot. In the 
remaining slots, there will be a certain rate loss due to the 
variations in the interference footprint as a result of other cells’ 
optimization. In order to minimize this loss, we will introduce 
a sum-power constraint in the optimization problem (detailed 
in the next section), such that the overall interference caused 
by one cell to others is upper bounded. Moreover, as it will 
appear, the solution to the optimization problem takes the form 
of a water filling algorithmQ It is then reasonable to expect 
that most of the power will be often assigned to users in the 
cell center, which reduces substantially the variations in the 
interference power as seen by other cells. These observations 

’The water filling algorithm is a general power assignment algorithm for 
multichannel systems with a sum-power constraint, which will assign higher 
powers to channels with better conditions. 











IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. X, NO. XX, XX 201X 


4 


are verified by the simulation results in Section [V] which 
shows that this scheduling mechanism can achieve nearly the 
same performance as the joint (optimal) power allocation after 
one scheduling round, but with a much lower complexity. 
This whole mechanism consists of two procedures: planning 
the scheduling sequence, and designing an optimal power 
allocation scheme for each cell. The first problem can be 
solved by frequency reuse patterns, and the second problem 
is our focus in this paper. 

III. Uplink Power Allocation Schemes 
Based on the per-cell optimization approach illustrated 
above, in this section, we aim to find an optimal power 
allocation scheme to maximize the uplink sum rate in the 
central cell, which can be obtained by solving the following 
optimization problem^ 

N 

PnO = arg ^ max (14) 

c=l 

* 

where P“q denotes the optimal p“q, and the power allocated 
to each user is adjusted at the time scale of slow fading. 
Unfortunately, it is difficult to derive exact solutions for this 
optimization problem with the expectation in (|7]i. We note that 
a closed-form expression for the achievable uplink rate with 
ZF receivers has been given in [26], however, this expression 
is rather involved and does not readily facilitate the design 
of an optimized power allocation policy. Instead, we derive 
tractable lower and upper bounds for the achievable uplink 
rate, as well as an approximation which lies between these 
two bounds. Based on these three expressions, we then present 
their corresponding power allocation strategies obtained as 
solutions to ([T4li . which are simple and easy to compute. First, 
we consider the lower bound. 


A. Lower Bound 


1) Closed-form Lower Bound Expression: 

Theorem 1: The achievable uplink rate of the nth user in 
the central cell is lower bounded as 


pul \ pul,L 
-^nO — -^nO 


l0g2 


1 + 


Kl)/3ono(M - N) 


\ 


ixlH 


(15) 


L N 

E E pfiPoci + 1 

;=ic=i / 

Proof: See Appendix U □ 

This lower bound is much simpler than the original uplink 
rate and facilitates the design of the following power allocation 
strategy. 

2 ) Power Allocation based on Lower Bound: With (115b . the 
optimization problem ([l4li can be simplified by replacing the 
exact uplink rate with the lower bound as follows 


N 


Pm = arg 




(16) 


c=l 


^As stated in Section IlLCl a sum-power constraint (rather than per-user) is 
imposed to control the interference level. 


The following theorem gives the solution to this problem. 

Theorem 2: The solution to the uplink power allocation 
problem (ITbl l i^ 

where is chosen to satisfy E^^iPnO = and d„o — 

/?o„o(M - 7V)/(Ef=i Ef=iF“[/3oci + l). 

Proof: To maximize the sum rate in (ITbl) . it is obvious 
that the total power needs to be set to the largest value 

and each user’s power should be nonnegative. Using the 
Lagrange multiplier approach associated with Karush-Kuhn- 
Tucker (KKT) conditions, we can obtain the solution to (ITbl) . 
which yields (ITtI i. □ 

Note that this power allocation is related to the number of 
BS antennas, the number of users and the large-scale fading 
coefficients. In particular, as M oo, 1/dno —t 0, and the dif¬ 
ferences in the allocated power among users vanish, meaning 
that this allocation strategy tends to an equal power assignment 
with asymptotically large M. A similar conclusion given in 
[27] shows that as the signal-to-noise ratio (SNR) grows, the 
power allocation that maximizes the sum rate converges to an 
equal power assignment. In our case, this result is obtained by 
keeping the total SNR fixed, but increasing the number of BS 
antennas. These observations indicate the equivalence between 
the two asymptotic regimes, large number of BS antennas and 
high SNR, which is consistent with the power-scaling results 
given in [15]. 


B. Upper Bound 

1) Closed-form Upper Bound Expression: 

Theorem 3: The achievable uplink rate of the nth user in 
the central cell is upper bounded as 


RnO < Ko“ = l0g2 ( 1 + PnOpoMM - N + 1) 
Qi^oyrhi^o) ( 


X E E V4(a) 

h=i i 

/1 \ 

e 

\ n/u \ / 

m=0 


U - 1)! 




where 


^OO g-xt 

El (x) = / - dt, Re(a;) > 0, 

Ji t 


(18) 

(19) 


is the exponential integral function of order 1 , Cofe — 
pIIPocI {I = 1,...,L;c= with k = N{1 - 1) -|- c. 

Moreover, Ao = diag (Coi, C02, • ■ ■, Cot), while T = NL, 
glyAa) is the number of the distinct diagonal elements of .Ao, 
Co(i> > Co(2> > • • • > Co(e(Ao)) are the distinct diagonal 
elements in decreasing order, Th{Ao) is the multiplicity of 
Co{;i), and \hj{Ao) is the (/i, j)th characteristic coefficient of 

A [ 28 ], 


^Observe the water-filling form of the solution. On the uplink, users get their 
transmit power information through the feedback from the BS, which acquires 
the channel information (lai'ge-scale fading) for all users and determines the 
uplink transmit powers. Here, the user scheduling involving fairness is not 
considered, but will be considered for future work. 
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Proof: See Appendix HIl □ 

We can now use this upper bound to provide a correspond¬ 
ing power allocation strategy. 

2) Power Allocation based on Upper Bound: After replac¬ 
ing the exact uplink rate R'^q in (fl4li with the upper bound 
given in (fT^ . the central cell optimization problem 
becomes 

N 

Ko = arg ^ max ^ ( 20 ) 

c=l 

and the solution is presented in the following theorem. 

Theorem 4: The solution to the uplink power allocation 
problem (l 20 l i is 

where /ig^ is chosen to satisfy '^^=1 Pnb — and 


Qi-^Ophi^o) ( 




fcno = ^o™o(M - A + 1) ^ Y.\ Vj(-4o) 

h=l J=1 I \J )■ 

j-2 


e«o<i.> El 






0(h) 


. ( 22 ) 




Proof: Follows the same steps as the proof of Theorem 


Observe that, despite the involved structure of (l22l i. the co¬ 
efficients kno (for different users) differ only through the factor 
/?Ono and the remaining factors can be regarded as constant. 
Moreover, as M —)> oo, 1 /kna —t 0, and the differences among 
users in this power allocation scheme disappear. Again, a 
tendency towards an equal assignment policy is observed with 
M growing without limit, which agrees with the observation 
in Theorem^ 

Next, we give a third tractable expression for the achievable 
uplink rate—a new approximation, which is proven to lie 
between the upper and lower bounds derived above. As shown 
later in this paper, our simulations demonstrate that this 
result is particularly accurate over a wide range of operating 
conditions. 


C. Approximation 

1) Closed-form Approximation: To obtain the approxima¬ 
tion, a key useful tool is given first. To the best of our 
knowledge, this approximation is new and provides a useful 
general tool for studying ergodic capacity. 

Lemma 1: If X and Y are independent positive random 
variables, then 


E 




log 2 


E{y}; 


(23) 


where 


logs 



< logs 


^{y}) 


< log2 




(24) 

Proof: See Appendix Hill □ 

Note that the approximation in this lemma lies between 
two bounds obtained by a direct application of Jensen’s 


inequality. Therefore, we can conclude that log 2 ^ 

is a relatively tight approximation of E{log 2 (l + a)} 
least tighter than one of the two bounds). It is also important 
to note that the lower and upper bounds given respectively in 
Theorems Q] and |i] are obtained by applying the correspond¬ 
ing bounds in (l24l i. Moreover, this approximation becomes 
asymptotically exact as shown in the following corollary. 
Corollary 1: Let X and Y in Lemma Q] be given as 

X = 1/Y^'\x,,Y = (25) 

where Xi (i = 1,..., Of) and yj (j = 1,..., Of) are both se¬ 
quences of positive square-integrable random variables which 
are not necessarily independent across i and j, respectively, 
and as 6*1 oo, —>■ u; G (l,oo). Moreover, there exist 

ei > 0 and £2 > 0 , such that 


liminf E {xi} > ei, liminf E {i/j} > £ 2 - (26) 

i j 


Then, as 0i 00 , 


E{log,(l + J)}-log,(l + |gl)^ 0 . ( 27 ) 

Proof: See Appendix IIVI □ 

When Lemma [7] is applied to the achievable rate analysis, 
9i denotes the number of BS antennas. Therefore, according 
to Corollary [7] the approximation (|23]) will be tighter as the 
number of BS antennas increases and becomes asymptotically 
exact as M —J- 00 . Note that a similar corollary can be easily 
obtained to show the asymptotic exactness of (l23l l for the 
case X = which becomes useful when analyzing 

maximum ratio combining receivers. 

Based on these results, the following theorem gives a 
tractable expression for the achievable uplink rate in (|7]i. 

Theorem 5: The achievable uplink rate for the nth user in 
the central cell is approximated as 


/ 


\ 


R. 


nO 


Ko = log2 


1 + 


ptoPono{M-N+l) 


(28) 


L N 

E E pfiPoci + 1 

\ /=lc=l / 

Proof: Divide the numerator and denominator of (ITJi by 
||ao„o||^- Then, 


^nO — IE \ log2 


1 + 


PnoPonO ZOnO 

E E Pd l^odl^ +1 

1=1 C=1 


\] 


)) 


(29) 


From the proof of Theorem [7] we know that zo„g and ^g^ are 
independent. Hence, applying Lemma [7]in ( |29] | yields 


^ 6 ui 


*^n 0 


Ko = l 0 g 2 


1 + 


Ln0/30n0E{z0n0} 


E Ep:;[e{icoc/|H 

1=1 c=l '■ 


+ 1 


(30) 


Plugging drol l and (l46l) into (l30l l leads to the desired result. □ 
Note from Lemma |7] that the approximation in dTHI ) lies 
between the lower and upper bounds in (ffSl) and dTSl l. and as 
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shown later in the numerical discussion, this approximation 
is more accurate than the aforementioned bounds. Moreover, 
according to Corollary\I\ this approximation becomes progres¬ 
sively more accurate as the number of BS antennas increases. 
Next, we will use it for the power allocation analysis, 

2) Power Allocation based on Approximation: With the 
tractable approximation (|28] |. the central cell optimization 
problem becomes 

N 

PnO = arg ^max '^Rco, (31) 

c=l 

and the following theorem presents the solution. 

Theorem 6: The solution to the uplink power allocation 
problem (l3Tl i is 

(32) 

where satisfies X]n=iKo = and — 

/?o„o(M -N + l)/(Eti + l). 

Proof: Follows the same steps as the proof of Theorem 

m □ 

When M ^ N, the approximation in (l28l l approaches the 
lower bound in (flSl l asymptotically (as M — N oo) . 
Therefore, for M ^ N, the power allocation scheme based 
on the approximation is also equivalent to that based on the 
lower bound. Hence, the power allocation strategy in Theorem 
|6] also tends to an equal power assignment as M ^ oo. 

We have provided three different power allocation schemes 
based on three tractable expressions for the achievable rate; 
an upper bound, a lower bound, and an approximation which 
lies between these two bounds. In Section |V] we will evaluate 
the performance of these three schemes with some numerical 
comparisons under different operating conditions, 

IV, Downlink Power Allocation Scheme 

Similarly with the uplink case, we assume the total power 
transmitted by one BS is constrained as 
Then, we aim to find an optimal power assignment scheme 
to maximize the cell sum rate. That is, we need to solve the 
optimization problem: 

N 

PnO = arg ^ max (33) 

c=l 

where P'^q denotes the optimal p^g, and the power allocated to 
each user is adjusted at the time scale of slow fading. The exact 
solution of this optimal problem is hard to obtain. Therefore, 
aiming to obtain a simple power allocation scheme, we first 
present a new lower bound for the achievable downlink rate. 


A. Lower Bound 

We first require the following preliminary result. 


Lemma 2: To satisfy the constraint (|9]l, the constant scalar 
ai for the ZF precoding matrix is given by 


ai = 


M-N 
Sn=l ^1 Plnl 


(34) 


Proof: See Appendix [V] □ 

We now give a tractable lower bound for the achievable 
downlink rate. 

Theorem 7: The achievable downlink rate of the nth user 
in the central cell is lower bounded by 


R. 


nO 


> <g^ = 10g2 


^ P^},(M-iV)/Agg 

V 1=1 PlclMl j 


(35) 


where Au = YJk=i l/Afcc 

Proof: See Appendix IVII □ 

Note that this lower bound is more tractable than the exact 
downlink rate (fOT l. Since the ZF precoding vector hid in (fOl) 
is different across cells and users, the approach to obtain the 
upper bound and approximation in the uplink case does not 
apply here. Therefore, in the following subsection, we use this 
lower bound to conduct the power allocation analysis. 


B. Power Allocation 

With (1351) . the optimization problem (l3^ can be simplified 
by replacing the exact uplink rate Rfk with the lower bound 

ndl.L 

PcO as 

N 

Pno = arg ^ max ^ (36) 

E pfo<P’'^ ‘==1 

c=l 

and the solution is given by the following theorem. 

Theorem 8: The solution to the downlink power allocation 
problem (1^ i^ 

Fno = (mo" - —) . (37) 

V S„o/ 

where satisfies Yln=iP^o = P'^^ and s„o = 

((M - N)/Aoo)/(Eti Eli Sef + l)' 

Proof: The derivation follows the same procedure as that 
shown in Theorem^ □ 

Note that the effective power allocation for the downlink is 
related to the number of BS antennas, the number of users and 
the large-scale fading coefficient. Similarly with the uplink, 
as M —>■ oo, 1/snO —5“ 0. Hence, the allocated power 
will approach P'^^/N, that is, the equal power assignment. 
Therefore, an important insight we can draw from comparing 
Theorem |?] and |6] is that regardless of uplink or downlink, 
the equal power allocation policy tends to be optimal as the 
number of BS antennas grows without limit. That means that 
large antenna arrays produce the same effect as increasing the 
transmit power and thus they can be used to cut down the 
transmit power for both uplink and downlink. 

“^The channel knowledge needed in the downlink power allocation (large- 
scale fading coefficient) is acquired from the uplink due to channel reciprocity. 
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V. Numerical Results 

In this section, we provide numerical results for a set of 
19 cells with radius Tc = 1000 meters in Fig. [T] N users 
are distributed randomly and uniformly in each cell, with the 
exclusion of a central disk of radius Vh = 100 meters around 
the BSs. The large scale fading is modeled using /3ici = 
ZiciliriciIrhY, where Zid is a log-normal random variable 
with standard deviation tr, v is the path loss exponent, and 
riel is the distance between the cth user in the hh cell and the 
7 th BS. In our simulations, we choose a = 8dB and v = 3.8. 
The small-scale fading is assumed to be Rayleigh distributed. 
Here, we consider a simple case where only adjacent cells 
produce interference (while further cells are neglected). Then, 
according to the scheduling approach described in Section 
Ill-CI the whole set of cells can be divided into three groups, 
marked as 1, 2 and 3 in Fig. [T] and each group conducts 
power allocation in a single time slot. Since this scheduling 
approach is based on the per-cell optimized powers, we will 
first focus on the performance of a single target cell (central 
cell in Fig. 1), during the first time slot in which interference 
conditions (i.e., the powers in the other interfering cells labeled 
2 and 3) are fixed. The network performance, taking into 
account dynamical effects of the interference footprint, will be 
considered subsequently, where we evaluate the sum rate of 
the entire 19 cells and compare our scheduled power allocation 
with the joint (optimal) power allocation. 


\ / 
% / 



Fig. 1: Simulation model. 


A. One-cell Performance 

1) Uplink: For the uplink, we first evaluate the accuracy of 
the three expressions for the achievable rate and the perfor¬ 
mance of the power allocation schemes based on them. Then, 
the best power allocation strategy is selected to quantify the 
benefits over equal power assignment in a range of scenarios 
with variable M/N ratios. As we will see, the approximation 
is the most accurate expression and the corresponding power 
allocation scheme can gain noticeable improvements over 
the equal power assignment for practical M/N ratios. The 
transmit powers of users in surrounding 6 cells are set to lOdB 
for the uplink analysis. 

In Fig. 13 the simulated uplink sum rate in 0 is compared 
with its lower bound in ( fTSl ). upper bound in ( fTSl ) and 
approximation in ( |28] |. Results are given under two different 
total transmit powers—20 and 30dB, and the transmit power of 


each user is assumed to be the same. Clearly, in all cases, the 
uplink rate grows with the increasing number of BS antennas. 
Moreover, we see a close agreement between the simulation 
results and our analytical lower bound and approximation. 
In particular, the approximation is almost indistinguishable 
with the simulated values and lies between the lower and 
upper bound, as expected from Section IIII-CI Observe that 
the lower bound and the approximation are very close to each 
other; indeed, this can be anticipated from (fTSl) and (l28l) when 
M » 



Fig. 2: Uplink sum rate per cell vs. the number of BS antennas with equal power 
allocation, where N = 10 users. 

To show the distinction between the lower bound and 
approximation, we present simulation curves for M closer to 
N in Fig. [3 Here, the simulated uplink sum rate as well as 
the bounds and approximation are plotted against the uplink 
total transmit power. We can see that, with M closer to N, 
the gap between the lower bound and approximation becomes 
more evident than that in Fig.|2] and in this case, the accuracy 
of the upper bound is comparable to that of the lower bound. 
Importantly, the approximation is still the most accurate. Next, 
we investigate the proposed power allocation schemes. 



Fig. 3: Uplink sum rate per cell vs. the total transmit power P'^^, with equal power 
allocation. 

Fig.|4]plots the different power allocation schemes described 
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in Section HIJ under both multicell and single-cell scenarios]! 
For multicell systems, the power allocation (fTTI i based on the 
lower bound, (l2Tl i based on the upper bound and ( |3^ based on 
the approximation are compared with the equal power assign¬ 
ment. Clearly, all power allocation schemes achieve noticeable 
improvements on the uplink rate, especially when the number 
of BS antennas is within a few hundred. As expected, these 
improvements become less significant as the number of BS 
antennas increases, and we observe that, when M ^ oo, 
the proposed power allocation schemes approach equal power 
assignments. It can be observed that the power allocation based 
on the approximation leads to the maximum enhancement, 
with a very slight improvement over the assignment based on 
the lower bound. This is a direct consequence of the increased 
accuracy of the approximation over the bounds. 



Fig. 4: Uplink sum rate vs. the number of BS antennas with different power allocation 
schemes, where — 20dB and = 10 users. 

To illustrate the effect of these power allocation schemes 
further, we define the relative gain of the power allocation as 
r] = (CpA — Ceci) /Ceq, where Cpa denotes the sum rate with 
our power allocation schemes, and Ceq denotes the sum rate 
under equal power allocation. 

Fig. El shows that, for the proposed power allocation 
schemes, significant gains are attained even with a few hundred 
antennas. For example, when M = 100, the power allocation 
based on the approximation in multicell systems can obtain 
almost 14% gain. Moreover, the gap in the multicell scenario is 
enhanced compared with that in the single-cell scenario. This 
is due to the addition of interference from other cells, which 
makes each user’s receive SNR smaller, thus leading to a better 
performance of the water filling algorithm. We thus anticipate 
this difference to grow as the number of cells increases. 

The above results indicate that the gain brought by the 
power allocation schemes diminishes as the number of an¬ 
tennas increases while keeping the number of users fixed, 
i.e., as M/N oo, which is sometimes referred to as 
the “massive MIMO regime”. However, this might not be 
representative of practical scenarios where (i) the number of 

^The power allocation schemes in single-cell scenaiio are obtained by 
setting L = 0. Then, the power allocation based on the approximation and 
that based on the upper bound have the same expression. 



Fig. 5: Relative gains of different uplink power allocation schemes, where — 20dB 
and N = 10. 

antennas is limited or (ii) the number of users is also large 
and is comparable with the number of antennas. In order to 
explore the effects of the proposed schemes in such scenarios, 
we next focus on the regime where M and N scale together 
with a fixed M/N ratio. Since the power allocation based on 
the analytical approximation (l32l i has been shown to attain 
the maximum gain, we only consider this strategy in the 
subsequent simulations. 

Fig. |6] shows the effect of optimized power allocation with 
different M/N ratios. We find that under a fixed M/N, the 
gain brought by the power allocation grows as the number 
of BS antennas increases, which differs from the conclusion 
under a fixed N. Since M and N scale together in this case, 
a larger number of antennas implies a larger number of users 
and, consequently, an increased interference which makes the 
receive SNR smaller. Notably, it is precisely under low SNR 
conditions where the power allocation optimization brings 
most benefit. The same arguments explain the fact that larger 
gains are attained at smaller M/N. It is also interesting to 
note that the curves in Fig. |6] are all approximately linear with 
no offset, which suggests that the relative rate gap between 
schemes with and without power allocation remains constant 
no matter the value of M. Therefore, the relative gain for each 
M/N ratio can be obtained from the slopes of these curves 
via linear fitting. These results are shown in Fig. |7] 

Fig.|7]presents how the relative gains of the power allocation 
change with the M/N ratio. We give results for four different 
total transmit powers—10, 15, 20 and 25dB. As expected, 
the relative gains reduce with increasing M/N or P'^^, which 
is consistent with the property of the water filling algorithm. 
Importantly, Fig. [T] allows determining the M/N ratios which 
bring gains over a particular (required) threshold, and these 
results can be used to decide whether the power allocation is 
justified. For example, if we set the tolerable minimum relative 
gain as 10%, the M/N ratio with transmit power 10 and 15dB 
can be high; for example, more than 20. This implies that 
even if there are more than 20 BS antennas per user, one still 
obtains more than 10% relative gain by using optimized power 
allocation. However, as the transmit power is boosted, the 
relative gains of power allocation will diminish. For example. 
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Fig. 6: Uplink sum rate vs. the number of BS antennas with constant M/N. 

for a transmit power of 20dB, relative gains exceeding 10% 
are only observed for 12 BS antennas or less per user (i.e., 
M/N < 12); whilst for a transmit power 25dB, these gains are 
achieved for only 4 BS antennas or less per user. For system 
configurations in which the relative gain does not exceed a 
prescribed threshold (e.g., 10%), one may favor equal power 
allocation, as a consequence of its minimal complexity. This 
is further explored in Table [H which indicates the maximum 
number of BS antennas per user (i.e., max value of M/N) 
required for two prescribed relative gain thresholds (10% and 
20%), for different transmit power constraints. These results 
demonstrate that, even when there are 10 times more BS 
antennas than users, and even for moderate transmit powers 
(e.g., 15dB), one can still obtain quite substantial gains (e.g., 
more than 20%) by performing power allocation optimization. 



Fig. 7: Relative gains of the power allocation based on approximation. 


TABLE I: Maximum M/N ratios with different uplink relative gain thresholds 
and total transmit powers 


Minimum Relative Gain 

^ul 

lOdB 

1.5dB 

20dB 

25dB 

10% 

93 

30 

12 

4 

20% 

29 

10 

4 

2 


2) Downlink: For the downlink, we first show the per¬ 
formance of the lower bound for the achievable rate. Then, 
the corresponding power allocation scheme is evaluated for 
two different kinds of users—central users and edge users. 
We find that although the power allocation strategy still tends 
to an equal power assignment when M/N oo, it obtains 
considerable gains with practical M/N values. The transmit 
powers in surrounding 6 cells are set to 30dB for the downlink 
analysis. 

In Fig. [8l the simulated downlink sum rate in (fTSl) is 
compared with the analytical lower bound in (l35l l under two 
different total transmit powers—40 and 50dB. The transmit 
power intended for each user is assumed to be the same. 
Clearly, in all cases, the downlink sum rates grow as the 
number of BS antennas increases and we can see a close 
agreement between the simulation results and our analytical 
lower bound. 



Fig. 8: Downlink sum rate per cell vs. the number of BS antennas with equal power 
allocation, where N = 10 users. 

Fig. |9] shows the effect of downlink power allocation based 
on the lower bound dJTl i. Results are given for central and 
edge users, which are separated by O.Src. The power allocation 
scheme has a noticeable rate improvement for edge users, but 
for central users the improvement is almost negligible. This 
is because the SNR of edge users is much smaller than that 
of central users and, as discussed above, the gains due to 
power allocation are most prominent at low SNRs. In addition, 
when the number of BS antennas increases, the improvement 
brought by power allocation also reduces, and as M —cx), this 
power allocation scheme tends to an equal power allocation. 
These observations are all in accordance with our analytical 
results. 

To further illustrate the improvement due to power alloca¬ 
tion for the downlink. Fig. [10] shows the relative gains of the 
proposed power allocation scheme over the equal power as¬ 
signment. This power allocation policy provides the maximum 
enhancement for edge users and, as anticipated, these gains 
are almost negligible for central users. As expected, the gains 
diminish with M, vanishing asymptotically as M ^ oo, which 
means that the equal power allocation tends to be optimal. 
Note however that the gains for edge users are still noticeable 
for antenna numbers up to a few hundred. For example, taking 
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Number of BS Antennas (M) 


Fig. 9: Downlink sum rate of central and edge users vs. the number of BS antennas with 
the proposed power allocation scheme, where — 40dB and = 10 users. 



Fig. 10; Relative gains of the downlink power allocation scheme, where — 40dB 
and = 10 users. 


M = 100, power allocation yields a 10% gain on the downlink 
rate for edge users, in contrast to a 2% gain for central users. 

Similarly to the uplink analysis, for the downlink, we are 
also interested in further exploring when the power allocation 
scheme is justified. Due to the precoding matrix constraint, the 
downlink sum rate does not show the linear behavior shown in 
Fig. |6]for the uplink (meaning that the sum rate depends on M, 
and not only on the ratio M/N). Therefore, we prescribe some 
specific values of M and N to perform this analysis, where 
we only consider edge users. Table HU a presents the maximum 
numbers of antennas for a minimum acceptable relative gain 
with different numbers of users and transmit powers. We can 
see that, as expected, when the total transmit power or the 
relative gain threshold increase, the maximum M decreases. 
Moreover, the maximum M grows as N increases and this 
is because more users result in more interference and lower 
SNR, which ultimately leads to a higher relative gain from 
the power allocation. In Table |II]b, we present the minimum 
numbers of users for an acceptable relative gain with different 
M and transmit powers. These two tables can be used to 
determine whether the downlink power allocation is justified. 
For example, when N = 10 with transmit power 35dB, the 
power allocation guarantees a minimum relative gain of 10% 


for any M up to 288. But for 40 and 45dB, such gains are 
achieved for M up to 99 and 38, respectively. Similarly, if M 
is set to 100, according to Table Hllb. a system with transmit 
power 35dB can obtain at least 10% relative gain with more 
than 7 users. But for 40 and 45dB, to achieve this target, 
there must be more than 10 and 16 users, respectively. As for 
systems that can not reach the prescribed minimum relative 
gain, the equal power allocation may be more suitable. 

TABLE 11: Relationship between M and N with different downlink relative 
gain thresholds and total transmit powers 



B. Multicell Performance with Scheduling Mechanism 

The previous results demonstrated the performance from 
the perspective of a given cell, assuming that the interference 
power originating from the surrounding cells was fixed. This 
allowed us to focus on the relative pros and cons of the 
different power allocation strategies under ranging conditions. 
Now, we adopt a network perspective, and evaluate the per¬ 
formance of a set of cells with our scheduling mechanism for 
power allocation, taking the uplink as an example. We adopt 
the power allocation strategy based on the analytical approx¬ 
imation (l32]) . as it was shown in Section IV-A.ll to yield the 
best performance among the different schemes considered. We 
apply this strategy for each cell and evaluate the sum rate of 
the 19 cells in Fig. [ijl The initial transmit power for each user 
is lOdB. In time slot i (for i = 1,2, 3,...), only cells within 
group j = mod (i — 1,3) -f 1 adjust their powers. Thus, during 
each scheduling round (comprising 3 time slots), each cell will 
operate with optimized power for one slot. In the subsequent 
two time slots, it will incur some performance loss due to 
a changing interference footprint caused by power allocation 
being performed in the surrounding cells. As a benchmark, 
we compare with the optimal performance achieved via joint 

®When evaluating the performance of these 19 cells, some adjacent cells 
out of this cluster, which bring interference to cells at the edge of this cluster, 
should also be considered. Since the number of power allocation objectives 
must be finite, we assume these out-of-cluster cells have constant transmit 
powers. 
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Fig. II: The uplink sum rate of 19 cells vs. the index of time slots, where M — 20, 
AT = 5 and = SOW. 


power allocation across all 19 cellsHland with the simple equal 
power allocation strategy (i.e., no power adaptation). 

From Fig. [TT] we find that, in slots 1 and 2, the joint 
method performs better than the scheduling one. But from 
slot 3, the gap between them is negligible. That is, after 
a single scheduling round, the scheduled power allocation 
achieves almost the same performance as the joint one, and this 
nearly optimal performance is maintained in the subsequent 
time slots. The joint method, highly complex and requiring 
high-level communication among BSs, entails a challenging 
implementation in practice. At a much lower complexity, the 
scheduled power allocation achieves nearly the same result 
(with differences inappreciable) after the first scheduling round 
which can be seen as a transient effect. Therefore, considering 
both complexity and performance, the scheduling mechanism 
is an attractive solution. 

VI. Conclusion 

We have characterized the gains in achievable rate brought 
by power allocation schemes in multicell MIMO systems with 
large antenna arrays. With ZF receivers, we have derived 
lower and upper bounds for the achievable uplink rate and 
an approximation which lies between these two bounds. With 
ZF precoders, we have derived a relatively simple lower bound 
for the achievable downlink rate. As opposed to a joint (across 
all cells) optimization, we have proposed a much simpler 
scheduling method to plan the power allocation arrangements 
for different cells. Based on this, new power allocation strate¬ 
gies have been presented, which can bring considerable gains 
over the equal power allocation assignment. These gains are 
substantial for moderate array sizes up to a few hundred, 
and for the case where N is allowed to scale with M, these 
gains increase as M grows with fixed M/N. Based on these 
results, we provide applicable values of M/N with a certain 
minimum relative gain in the achievable rates, which can be 
used as practical design rules to justify the use of optimized 

^The joint power allocation problem does not admit an analytical solution; 
thus, it was evaluated numerically using the function “fmincon” in the 
MATLAB optimization toolbox. 


power allocation schemes over the equal assignment policy. 
This reveals the applicability of our allocation schemes under 
a wide range of scenarios with practical numbers of users 
and antennas. From a network point of view, the proposed 
scheduling mechanism achieves almost the same performance 
as the joint power allocation after one scheduling round, with 
much reduced complexity. 


Appendix I 
Proof of Theorem[T| 

Application of Jensen’s inequality to (|7]i gives 


<0 > Ko" = log2 


/ 

1 +- 


PnO 


\ 


V 


E E PSE{|atooSoc;|^} +E{||ao„oin^ 

(38) 

which can be further written as 


/ 


Ko" = log2 


1 +- 


PnO 


\ 


, E E PSE{ICocifl|ao„of}-bE{||ao„oin. 

\ 1=1C=1 / 


(39) 

where ^oc/ = a^oSoci/llaonoll- Since god has a rotation- 

can be regarded as a column 


invariant distribution and -n- 

l|a_ 

of a rotation matrix, 4oc/ has the same distribution as the 
elements of gocz and is independent of ao„o- Therefore, 


1e{|^0cz|'} =/3ocZ, 
and ( |39] | can be expressed as 
( 


= log2 


1 + 


PnO 


(40) 


\ 


(^E EpSE{i?ocir}+ijE{iiao.oiny 


Ho^oHoo) 


(41) 

, which is chi-squared 


Let zono — 

distributed with probability density f29] 

p — ^OnO 

and therefore, E{l/zono} = 1/{M — N). Then, with 

||aOno||^ = 


= l//3on020nO, 


(Go^oGoo)' 

IE |||ao„o||^| = Ponoi.Nl — N). 

Plugging drol l and ( |4^ into (HTI i yields the result. 


(43) 


Appendix II 
Proof of Theorem[2] 

Applying Jensen’s inequality in (13l gives 


/ 


RnO < Ko = l0g2 


1 +E ^ 
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PnO 


\ 


E E PfMonOZOclt + llaonof , 

^ 1=1c=l > ' 

( 44 ) 
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Recalling the definition of and Zono in the proof of 
Theorem |7] and the fact that ^ocZ is independent of zono, we 
obtain 



/ 

f \ 

\ 

T^Ul.U 1 

Rn6 =log2 

1+E {pJJoPonOZOno} E < 

1 

L N 

> 


1 

E Ep“1Coc1 + i 

<l=lc=l > 

/ 


(45) 

Recalling (l42l i. we can get the following expectation 

E {p'^oPonOZOno} = Ko/3ono(M - iV + 1). (46) 

Now, we focus on the remaining expectation in ( |45] ). 
We know that pfi\^ocif ~ Exp(l/p“J/3oci), where Exp(a;) 
is the exponential distribution with parameter x. Let vo = 
pS ICocil^, which can be written as a single summa¬ 
tion, i,e,, Vo = V'ofc, where i/'ofc = Pci l^odf with k = 

N {1 -l) + c. Hence, V'ofc ~ Exp(l/Cofc), where Cofc = P^jPod 
with k = N{I — 1) + c. The probability density function of vq 
is given by [30, Theorem 2] 

e(-4o)Tfe(Ao) c-j 

/(«o)=E E (47) 

which leads to 
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By replacing vq with — s, 
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(48) 

s 

e‘’ 0 (h) ds. 

(49) 


With s-l 1 _ 1 = (s _ 1 ) ^ -f fil ^ + ■ • • + l), ( l49] l can be 

written as 


f 


3-1 -t2Q 
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‘ — oo 




fil ^ -f fil ^ -f • • • + 1) ds 


+ (-l)^ 

Using integration by parts, we can get 


1 r°° 1 

Jo s-l 
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m TTi+l 


and we also have 


f-OO ^ 

Jo s-l' 


e°^°‘ds = eil“Ei(l/a), 


Substituting (fSTl i and (l52l i into (fSOl l gives 

r°° , 4-1 -"0 . , 
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Jo uo-fl 
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The theorem follows by substituting (l5^ . (|48] | and 

dSll. 


(51) 


(52) 


(53) 


into 


Appendix III 
Proof of Lemma[T| 
From Jensen’s inequality. 


+ 0 + f) 

Since X and Y are independent. 
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Utilizing this on the left and right hand side of (l54l) . we can 
obtain (l24li . 

Comparing dH and (l24l) . it is obvious that 
log 2 lies between the lower and upper bound 

of E {log 2 (1 + y) }■ Therefore, we get the approximation 
expression in (l2Tt , 


Appendix IV 
Proof of Corollary[I] 
According to the law of large numbers, 

N ‘ T,. _ ^ \ 

Then, since ^ Dili Ejaii} is bounded away from 0, 
1 1 


- E ^ zii — ^ E ^ ® ^ 0, as 9i —>• 00 . (56) 
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Combining (|57] | and (ISST i, we get 
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Therefore, since ^E{F} is also bounded away from 0, 
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Vy “ E{y }) —^ 0 , as 01 ^ 00 . (63) 

By recalling the definition of almost sure convergence, we 
know that Ve > 0, there exists 60 , such that for 01 > 00, 
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where Pr (•) denotes probability. Hence, 
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and therefore 
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Then, Corollary [T] follows by a direct application of 
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Appendix V 
Proof of Lemma|2] 

Utilizing the property tr (BjiB^) = tr (B^B;;) on (| 9 ]l, we 


can obtain 

E{tr(B^Bii)} = 1. (67) 

Then, substituting (fTOl i into (l67l i yields afE |tr (G^G;*) ^| = 
1, which implies 

= yi/E{tr(G^G*)-i}. (68) 

We know that 

Ejtr (g^G„)“'| =tr (^E|(G^Gii)"'|^ (69) 

= ;^e|[(g^G„)“'1 |. (70) 

n = l L L JnnJ 

From ( |4^ . we have 

Applying this in (l6^ yields 

E(„(GSG„)-‘} = ^f (72) 

n=l 


Therefore, since tr (Gj^G*;) ^ = tr (Gj^Ga) o; is got by 
applying ( l72b in ( l68T l. 


Appendix VI 
Proof of Theorem[2] 
Applying Jensen’s inequality in (fTsT l yields 


/ 


7?n0 > Rn^ = l0g2 


1 + 


aoPno 


E EPSEllgFnoblclI^l 

1 = 1 c=l 


+ 1 


which can be further written as 


(73) 


R"r^<^ = logj 


/ 

1 + 


^ 2 dl 

^OPnO 


L N fl T u r 1 

EME +1 


. (74) 


From the proof of Theorem [7] we know that has 

the same distribution as go„o and is independent of hi^. 
Therefore 


= log2 


and 


1 + 


2 dl 


\ 


E + 


E||ginobici| /||bici|| |=/3(n0- 

With ||bicz||^ = (G^G*j) ^ and (fTTl i. we have 


(75) 

(76) 


lE{l|b(.iin 


a? 


pidiM-NY 


ill) 


The theorem then follows by substituting ( l34b , ( f76] l and (ITTT i 
into dTSl l. 
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